Anisotropic s- wave superconductivity in borocarbides LuNi 2 B 2 C and YNi 2 B 2 C 
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The symmetry of superconductivity in borocarbides 
LUN12B2C and YM2B2C is an outstanding issue. Here an 
anisotropic s- wave order parameter (or s+g model) is pro- 
posed for LuNi2B2C and YN12B2C . In spite of the dominant 
s- wave component the present superconducting order param- 
eter A(k) has nodes and gives rise to the \fH dependent spe- 
cific heat in the vortex state (the Volovik effect). This model 
predicts the fourfold symmetry both in the angular dependent 
thermal conductivity and in the excess Dingle temperature in 
the vortex state, which should be readily accessible experi- 
mentally. 

PACS numbers: 74.60.Ec, 74.25.Fy, 74.70.Dd 



I. INTRODUCTION 

The superconductivity in rare earth- (R) transtion 
metal borocarbides is of great interestM, in particular 
its interplay... with magnetism and superconductivity is 
fascinatingEiu. However in the following we limit our- 
selves to the nonmagnetic borocarbides LUM2B2C and 
YM2B2C . They have a relatively high superconduct- 
ing transition temperature T c =15.5 K and 16.5 K re- 
spectively. Although the dominance of the s-wave 
component in A(k) has been established by substi- 
tuting Ni by a small amount, of Pt and subsequent 
opening of the energy gapaEl, a number of pecu- 
liarities are not expected in a conventional s- wave 
superconductorO. For example the sfB. dependence of 
the specific heat in the vortex state indicates a super- 
conducting state with nodal excitations sjmaila*, to d- 
wave superconductivity in high T c cupratesuSEHIj. Fur- 
thermore the presence of de Haas van Alphen (dHvA) 
oscillations in the vortex state of LUNI2B2C clasra._to 
H=0.2H c2 suggests again nodal superconductivitjO'EjO. 
In a conventional s-wave superconductor. dHvA oscilla- 
tions would disappear for H<0.8H C 2EiEj. In addition 
the upper critical field determined for LuN^I^C and 
YNI2B2C for field direction within the a-b plane exhibits 
clear fourfold symmetry somewhat reminiscent to d- wave 
superconductorsollj, furthermore 1/Ti from NMR ex- 
periments shows T 3 pojwer law behaviour consistent with 
nodal superconductorst3. These experiments clearly in- 
dicate that A(k) in borocarbides has to be an anisotropic 
s- wave order parameter. Furthermore a) A(k) has to 
have a nodal structure with the quasiparticle density of 



states (DOS) N(E)~ \E\ for \E\/A <1, which gives 
\[H dependence in the specific heat of the vortex stateB 
b) the nodal structure has to have the fourfold symmetry 
within the a-b plane which is consistent with the tetrag- 
onal symmetry of the a-b plane. These two constraints 
appear to suggest almost uniquely 
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A(k) = -A(l + sin 4 tfcos(40)) 



(1) 



or s+g- wave superconductivity. Here <f> are the polar 
and azimuthal angles in k- space respectively. We show 
in Fig. 1 A(k) which exhibits clear fourfold symmetry. 
The four second order nodal points of A(k) are given 
by (i?,<£) = (§,±f) and (f,±x) which dominate the 
quasiparticle DOS at low energies: 



N(E) _n\E\ E 2 
N 4 A ( A> 



(2) 



where No is the normal state DOS. In constructing 
A(k) p-jpej-have made use of a similar approach as in 
MgB 2 BEEI In the s+g model of Eq.(g) we assume 
the equality of s and g amplitudes to have N(E)~ \E\ 
down to lowest, energies. Recent thermal conductivity 
measuremental3 report a gap anisotropy of at least a 
factor of 10, the fine tuning of s and g amplitudes in 
Eq.(|l]) therefore has a tolerance of 10%. There is no sym- 
metry reason why the amplitudes (or pair potentials) of 
inequivalent representations like s and g should be very, 
close. However from the bandstructure of borocarbideslia 
it may be argued that the pair potential at the nodal 
points given above is indeed strongly suppressed. The 
main Fermi surface sheet shows lobe like structures along 
the [110] directions which have strong nesting with a wave 
vector parallel to a. This wave vector appears as the in- 
commensurate ordering vector in the magnetic borocar- 
bides (Lu,Y replaced by rare earth elements). Therefore 
the lobe states at (■&, (£)=(§ , ±f ) and (§,±^) tend to 
an instability in the particle-hole channel which strongly 
depresses the effective potential and associated A(k) for 
Cooper pairing at these points. The approximate fine 
tuning (up to 10%) of s and g amplitudes may be caused 
by this peculiar Fermi surface feature of the borocarbides. 

In the following we shall first consider thermodynam- 
ics and transport of the borocarbides for zero field for 
the proposed gap function. Then we will study the field 
angle dependence of specific heat and thermal conductiv- 
ity which exhibit the fourfold symmetry. We apply the 
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FIG. 1. Normalized gap function f(k)= A (k) /Aofthes+g 
model. 



same technique developed in Ref. [ [7| |22j, |23|, |24j . Also we 
predict the fourfold symmetry in the excess Dingle tem- 
perature in the vortex state in borocarbides in a planar 
magnetic field. 



II. THERMODYNAMICS AND TRANSPORT 
PROPERTIES 

First of all A(k) =Af(k) given in Eq.(0) leads to the 
quasiparticle density of states 

"W- 1 fdOBe - N 
JV 4tt J ^/x 2 - f 2 



\x\ f dyF(y)Re^=^= 

Jo v x ~~ y 



(3) 



where x=E/A and 

F(y) = ± -= 



dz 



with uo = (1 - V|l-2j/|)3 (4) 

we note that F(l-y)=F(y) holds. The quasiparticle 
density of states is evaluated numerically and shown in 
Fig.2. For \E\/A <1 we obtain 



N 4 A 4tt A '"' 
then the low temperature specific heat is given by 



(5) 



FIG. 2. Quasiparticle density of states. Logarithmic sin- 
gularity occurs at E= A due to the saddle points at i? = 0, 7r. 
The cusp at E=A is due to the gap maxima at (#, <^)=(?,0), 
(f,±f)and(f,7r). 



1N T ~ 4/ W( A j+ 80 1 A j 



(0) 



where is the Sommerfeld constant. Similarly the 
spin susceptibility and the superfluid density are given 
l.v' 



2 A v ' 16 V A ; 



Xn 
Ps(T) 



X 



= 1 - 

Ps(P) XN 



(7) 



Likewise the electronic thermal conductivity of the s+g 
model at low temperature is obtained in a universal form 

as 



K 

T 



it 2 n 



8 mA ® 

The prefactor 7r 2 /8 is specific for the s+g model. Here 
n, m are the electronic density and mass respectively. 
This is equivalent to k/k„ =3r/8A where K n is the nor- 
mal state thermal conductivity and T the quasiparticle 
scattering pate. The linear T behaviour of k has recently 
been foundO in LUM2B2C from which we extract T/A < 
0.02. 



III. ANGULAR DEPENDENT SPECIFIC HEAT 
AND THERMAL CONDUCTIVITY 

We are proposing that the angular dependent spe- 
cific heat and especially thermal conductivity in the 



2 



m: 



e window to look for the 
Indeed from the latter 



vortex state provides 
symmetry of A(k) E. 
Izawa et al have succeded in deducing the symmetry 
of A(k) in S^RjaQJHS, CeCoIn 5 EZI and more recently k- 
(ET) 2 Cu(NCS) 2 E3. First of all we have to construct the 
equation for the residual density of states in the presence 
of impurity scatterings. 



g = Re( 



Co 



V(C -ix)2+f2 



+ xtan 1 ( — ) 



(9) 



where Co=lm\j^o -^?ri(u)/A) with uj giving the renor- 
malized frequency and x= |v- q|/A ~ | sin(0± j)\. Here 
2q is the sum of the pair momentum associated with 
a supercurrent circulating around each vortex and v • q 
is the Doppler shift connected with it. In the first line 
the brackets mean averaging over both Fermi surface 
and vortex lattice, in the second line the former is eval- 
uated up to the ± summation and the latter still re- 
mains. In the superclean limit defined by Co <C (%) or 
r < v a \[eH «T« A Eq. (ft) gives 



9 = j{x) 



2V2A 



1(6) 



(10) 



where v = y/v a v c and 1(9)— max (| sin#|, | cos 9\) for 0< 
6 < § • The function 1(9) is shown in Fig. 3. Here v a 
and v c are Fermi velocities in the a-b plane and along the 
c-axis respectively. The magnetic field is applied within 
the a-b plane at an angle 9 with respect to the a axis. 

In the clean limit with Co (x) or v a \feH <C T <C 
T <C A on the other hand we obtain 



9 = 9 



(0)( 
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vHeH) H A 



1 



32rA 



(l-cos(40))]) (11) 



From these expressions the field angular dependent 
specific heat in the vortex state may be derived. In the 
superclean limit we obtain 



InT V2A 



■1(9) 



(12) 



In the clean limit on the other hand the above 9 de- 
pendence is replaced by 



C„ 
7nT 



v>(eH) [H A 



32rA 



1 



(1 



cos(W))\ 
(13) 



where g(0)=N(0)/No in the absence of magnetic field. 

The thermal conductivity tensor in the vortex phase 
has been calculated inEj and in a planar magnetic field 
it is given by 



3 i 2 (eH) 



G4 



sin(20) 



(14) 



in the superclean limit and 



K 



rv X y 
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(15) 



i2rA 



in the clean limit. Here is i^xx 

(H=0). Therefore 

we expect the fourfold symmetry in the thermal conduc- 
tivity in the vortex state should be readily accessible in 
future experiments. On the other hand has recently 
been measured for field oriented along cO. In this case 
a similar calculation in the superclean limit for H^; H C 2 
leads to 



-.(H) = 3 v 2 a (eH) H-H cl 
c„ 7T A 2 (0) ~ H c2 (0) 



(16) 



pTjhis behaviour was indeed experimentally observed 
ir£j. In the clean limit k xx (H) is no longer exactly linear 
but has a logarithmic correction in H. Since T/A <0.02 
for LUM2B2C we can use the above equation for the su- 
perclean limit except for very small fields. 



IV. EXCESS DINGLE TEMPERATURE 

It is well known that dHvA oscillations can be seen in 
the vortex state as well when the quasipaipUcle|-d]amping 
is much less than the cyclotron frequencjO'E2lEil. How- 
ever in conventional s-wave superconductors the dHvA 
oscillation becomes invisible when H<0.8 H C 2. Therefore 
if dHvA oscillations. are seen even for H~0.2H C 2 as in the 
case of LuNi2B2CE3, this can be taken as a signature of 
a nodal superconductor. Since the excess Dingle temper- 
ature in the vortex state is due to quasiparticle damping 
caused by the Andreev scattering it should also exhibit 
the fourfold symmetry of the order parameter. The ex- 
cess damping due to Andreev scattering is evaluated as 




(17) 



where we defined 



J(d) 



1 
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dd(l +sin 4 i9cos(46»)) 2 
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= i (l + 3 cos(40) + — cos 2 (40)) 



(18) 
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FIG. 3. Angular dependence of specific heat C s ~ 1(0) and 
excess Dingle temperature ~ 3(0) in an external field in the 
a-b plane. is the angle between field direction and a- axis. 



That is we average A 2 (k) on the Fermi surface sliced 
perpendicular to H. The angular dependence of 3(9) is 
shown in Fig. 3 together with 1(9). In particular we 
find I(|)/I(Q)= -j= and J(f)/J(0)= 0.2587. The excess 

damping is reduced by a factor of \ for H || [1,1,0] as 
compared to the one for H || [1,0,0] for example. 



V. CONCLUDING REMARKS 

Here we propose a simple model for A(k) for nonmag- 
netic borocarbide superconductors with fourfold symme- 
try. The angular dependence of the specific heat, ther- 
mal conductivity and the excess Dingle temperature are 
worked out with this model. We hope that this work will 
stimulate further experiments on borocarbide supercon- 
ductors. 
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